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Abstract. We present measurements of exclusive tv + '° and rj production in pp reactions at 1.25 GeV and 
2.2 GeV beam kinetic energy in hadron and dielectron channels. In the case of 7r + and ty°, high-statistics 
invariant-mass and angular distributions are obtained within the HADES acceptance as well as acceptance 
corrected distributions, which are compared to a resonance model. The sensitivity of the data to the yield 
and production angular distribution of Zi(1232) and higher lying baryon resonances is shown, and an 
improved parameterization is proposed. The extracted cross sections are of special interest in the case of 
pp— ^pp?7 , since controversial data exist at 2.0 GeV; we find a = 0.142 ± 0.022 mb. Using the dielectron 
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channels, the 7r and rj Dalitz decay signals are reconstructed with yields fully consistent with the hadronic 
channels. The electron invariant masses and acceptance corrected helicity angle distributions are found in 
good agreement with model predictions. 

PACS. 13.75Cs 25.40Epl3.40Hq 

1 Introduction 

Meson production in nucleon-nucleon reactions in the kinetic beam energy range 1-2 GeV tests an important sector 
of strong interaction on the hadron level. It is the subject of extensive studies by both experiment and theory with 
the aim to establish eventually the link of hadron physics to QCD as fundamental theory. Most of the available high- 
statistics data originate from recent near-threshold measurements, at excess energies < 150 MeV, performed at the 
SATURNE, CELSIUS and COSY facilities. Phenomenological models, usually adjusted to data of these reactions, 
serve as step towards a concise theoretical description. Such models are based on, for instance, the one-boson exchange 
(OBE) approximation describing the production amplitudes by a coherent sum of meson, nucleon and baryon-resonance 
currents. They reveal that meson production is a complex process with important contributions from nucleon-nucleon 
final-state interactions and interferences between the contributing reaction channels even within a tree-level approach. 
This situation often leads to ambiguous model descriptions of the experimental results (for a review see pQ). For 
higher energies, i.e. excess energies > 150 MeV, the data base is more scarce and originates mainly from low statistics 
bubble-chamber experiments [2J. 

The production of n and rj mesons in nucleon-nucleon collisions is of particular importance because of their coupling 
to baryonic resonances. Hence, experimental data on one-7r and one-?? production provide quantitative information on 
hadronic interactions, as well as resonance excitations and resonance properties. 

Cross sections for pion production in the beam energy range from 0.6 to 1.5 GeV (excess energies between 140 
and 500 MeV) have been provided in the past by many experiments [3,4.5,6,7,8]. Meanwhile, also the precision of 
the measurements of differential distributions, which is essential for unravelling the reaction mechanism, has been 
improved [9, 10l [Tll[T2l[TB1ll41ll5] . The dominance of the intermediate Z\(1232) production and the peripheral character 
of the reaction, which increases with energy, stand out very clearly in this region. Comparison of the shape of various 
differential distributions from both exclusive reactions pp— >-pp7r° and pp— ^pn7r + in the range 0.6 — 0.94 GeV to 
calculations within the one-pion exchange (OPE) show a nice agreement [TUlIll) . However, the magnitude of the cross 
sections for both reactions, regardless of the choice of the form of 7r-nucleon interaction, is explained by the models 
within an accuracy of 20 - 30% only. This points to some missing elements in the assumed reaction mechanism, as 
for example the exchange of heavier mesons, contributions of heavier resonances or/and non-resonant ir production, 
as demonstrated by a recent detailed analysis of pp— >pn7r + and pp— >pp7r° reactions at a beam energy of 0.94 GeV 
within the framework of a partial- wave analysis |15) . 

At beam energies higher than 1.5 GeV, the interpretation of the low-statistics bubble-chamber data [TMTTlfTB] was 
based on the isobar model assuming an incoherent sum of contributions from the decays of various baryon resonances 
into pions. However, large uncertainties remain, due to the limited statistical significance of the corresponding experi- 
mental results. Further progress in the understanding of meson production in p+p interactions thus requires new high 
statistics data. 

For the one-r; production, precise data have been collected close to the reaction threshold (excess energies below 
120 MeV) Ppjpip2P5ptp5P5P7P5] and compared to various OBE models (for a survey cf. [I]). Most of the 
calculations indicate a dominant role of resonances, in particular the AT* (1535) formed via the exchange of virtual 
pseudoscalar (71", rj) and vector (p, uS) mesons. However, the models differ in the description of how the resonance is 
excited, which is due to the uncertainty in the nucleon-meson-V* couplings. Better constraints can be obtained from 
differential distributions and polarization observables, as demonstrated in [29]. The dominance of A"* (1535) seems also 
to persist at higher beam energies (2 — 3 GeV), as shown by a detailed analysis of the pprj Dalitz distributions by the 
DISTO collaboration [3D]. The latter work, however, does not provide absolute cross sections which are very important 
for the quantitative evaluation of the role of resonances. 

Under the assumption that intermediate baryon resonances play a dominant role in n, rj and p production, a model 
was developed [3Tj based on an incoherent sum of various resonance contributions. The matrix element of the Z\(1232) 
production was calculated within the OPE model [35], which had been adjusted to available differential distributions 
of pion production in the pp— >pn7r + channel at incident kinetic energies in the range 0.9-1.5 GeV. The other matrix 
elements were kept constant and were determined by fitting the total meson production cross sections. 

Meson production is an important ingredient of microscopic transport models which were developed to describe 
heavy-ion collisions. Such approaches rely on a realistic treatment of elementary nucleon-nucleon and meson-nucleon 
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interactions, as described in the previous paragraphs, to calculate double-differential cross sections of the produced 
particles [33l[34l[35][36][37l[38] . 

Electromagnetic decays of mesons and baryon resonances are sources of e + e~ (dielectron) pairs which play a 
prominent role in heavy-ion physics as penetrating probes of nuclear and hadronic media. Therefore, data on 7r, rj, and 
baryon resonance production in proton-proton interactions are essential for the understanding of e + e~ pair production 
in p+p, p+nucleus and nucleus- nucleus collisions. Recent precise measurements of the HADES collaboration underline 
the need to understand elementary sources of e + e~ pairs for the interpretation of heavy- ion data [39,40,41,42 . It turns 
out that the baryon resonances are especially important. They contribute to the dielectron spectra via direct Dalitz 
decay, A — >Ne + e _ , and via two-step processes in which the resonance decays into a nucleon and a meson with a 
subsequent e + e~ pair produced in the Dalitz decay of the meson (e.g. Z\(1232) ^Ntt° followed by tt° — > r ye + e~ or 
N(1535)— >N r\ followed by r\ — > 7e + e~) or two-body decays of the produced mesons (e.g. N(1520)— >-Np followed by 
p — >c + e~). 

At higher beam energies [43], the dielectron yield is a complicated cocktail resulting from decays of many mesons 
and baryon resonances, but the main dielectron sources in nucleon-nucleon collisions are still the tt° Dalitz decays 
for dielectron invariant masses M ee < 0.135 GeV/c 2 , the ij Dalitz decay contribution for 0.14 GeV/c 2 <M ee < 0.547 
GeV/c 2 and Dalitz decays of baryon resonances and light vector mesons (p and w) for M ee > 0.6 GeV/c 2 . 

In the past, experiments studying p+p and 7r+p interactions have either analyzed the hadronic or leptonic final 
states. The HADES apparatus [33] allows for the first time to measure hadron and e + e~ pair final states simultaneously 
with high statistics. With such data, it is possible to achieve a consistent description of meson production in p+p 
reactions in the hadron as well as in the dielectron channel. In this article, we present the first step in this direction 
and compare experimental results from the analysis of three reaction channels pp— >-pn7r + , pp— ^pp7r°, and pp— >-pp?7 
measured at kinetic beam energies of 2.2 GeV and of the corresponding dielectron channels pp — > pp7r° — > ppe + c~7 
and pp— — s-ppe + e - 7 obtained at a kinetic beam energy of 2.2 GeV. At 1.25 GeV, which is below the threshold of 
r\ production in pp reactions, only the reactions involving pions were analyzed. As reference model, we use predictions 
of the aforementioned resonance model of 31 , complemented with experimental results [3D]. Differential spectra have 
been measured with high statistics providing strong constraints on the production mechanisms as well as on the 
different resonance contributions. 

The polarization of virtual photons has never been measured in the Dalitz decays of the pseudoscalar mesons, 
although the distributions of the respective angles of the emitted lepton with respect to the virtual photon direction 
(the so called helicity) have been predicted for several sources [32] . It has also been suggested that such distributions can 
be used as important "fingerprints" to distinguish between different dielectron sources in the inclusive measurements. 
Indeed, as shown in (46) , virtual photon polarization appears to be a very important characteristics of the dielectron 
excess radiation originating from the hot and dense hadronic matter created in heavy-ion collisions at SPS energies. In 
contrast to the situation at the SPS, similar measurements at SIS18 energies show a sizable anisotropy in the helicity 
distributions [42] , Therefore, it is desirable to determine the relevant distributions for the ir° and r\ mesons, as well 
as for baryon resonances, such as the Z\(1232), and bremsstrahlung in p+p collisions. In this work, we present such a 
measurement for the r\ Dalitz decay which has been isolated in the pp— >ppe + e _ 7 reaction channel. 

Our paper has the following structure. The experimental set-up and event reconstruction will be briefly described 
in sec. [5] The data analysis methods and simulation tools are presented in Sees. |3Jand|31 Section [5] is devoted to the 
discussion of the results. We draw conclusions in sec. [5J 

2 Experimental aspects 
2.1 Detector overview 

The High Acceptance Dielectron Spectrometer (HADES) consists of six identical sectors covering polar angles between 
18° and 85° and between 65 and 90 % of the azimuthal range. While a detailed description of the set-up can be found 
in |44j . we summarize here only the features relevant for the present analyses. Proton beams with intensities up to 
10 7 particles/s were directed to a 5 cm long liquid-hydrogen target of 1% interaction probability. The momenta of 
the produced particles are deduced from the hits in the four drift chamber planes (two before and two after the 
magnetic field zone) using a Runge-Kutta algorithm |44j . The momentum resolution is 2-3% for protons and pions and 
1-2% for electrons, depending on momentum and angle [33! ■ The trigger for the HADES experiments consists of two 
stages: The first level trigger (LVL1) is built on different configurations of hit multiplicity measurements in two plastic 
scintillator walls for polar angles larger (TOF) and smaller (TOFINO) than 45°, respectively. The second level trigger 
(LVL2) selects candidates defined by a Ring Imaging Cerenkov detector (RICH) and information from TOF and 
an electromagnetic shower detector (Pre-Shower) behind TOFINO. The analysis of hadronic channels was based on 
LVL1 triggered events selected by either of the two following configurations: The first required a coincidence between 
two hits in two opposite sectors of the time-of-flight detectors, with at least one in the TOFINO. This configuration 
was optimized for the selection of pp elastic scattering events and also used for the pp— >pp7r° and pp— >pnir + reactions. 
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Fig. 1. (Color on-line) Correlation between momentum (p) and velocity (v/c) for particles without a RICH signal and with 
curvature corresponding to a positive charge. The overlaid yellow curves indicate graphical cuts to select protons and pions. 

The second configuration was based on a charged particle multiplicity of four or more, at least one hit in the TOFINO 
and at least two signals in opposite sectors. This selection was used for the reconstruction of pp— >-pp?7 hadronic and 
dielectron channels (i] — > 7r + 7r~7r° and r\ — > 7e + e~) as well as pp— >pp7r° dielectron channels (tt° — > 7e + e~) in 2.2 GeV 
collisions. For the 7T° and ij electromagnetic Dalitz decay measurements, in addition to the latter LVL1 configuration, 
a LVL2 decision was requested, i.e. at least one electron candidate in the RICH. 

2.2 Event reconstruction 

Five different final states were used for the study of exclusive 7r and r\ production in pp reactions, as summarized in 
table [TJ An important feature of the HADES apparatus, which is exploited in the present analysis, is the ability to 
measure both hadrons [4"Tll4"8] and electrons [39,40,41] in the same experimental run. Of special importance for the 
analysis of the dielectron channels is the suppression of tracks produced by photon conversion and consequently the 
reduction of the combinatorial background [44]. This is achieved using criteria related to the track quality and the 
distance and opening angle between neighbouring tracks. Finally, only e + c~ pairs with an opening angle larger than 
9° were propagated to the physics analysis. The remaining combinatorial background is subtracted from the measured 
unlike sign pair yields, using the arithmetical mean of the like sign pair (e + e + and e~e _ ) yield in the same event. 



reaction 


meson decay 
channel 


measured 
exit channel 


incident 
energy 


pp— !>pn7r + 




7r+p(n) 


1.25, 2.2 GeV 


pp^ppvr 


all 


PP (*") 


1.25, 2.2 GeV 


7r u — > 7e + e~ 
(BR 1.12%) 


ppe + e (7 ) 


2.2 GeV 


pp^ppjy 


Tj — > 7T + 7r~7T U 

(BR 22.7%) 


PP7T + 7T~ (7T U ) 


2.2 GeV 


77 — > "/e + e~ 
(BR 0.68%) 


ppe + e" (7) 


2.2 GeV 



Table 1. Investigated reaction channels for exclusive it and r\ production reactions. 



Since the RICH is hadron blind in the given energy range and the particle multiplicity is low, electrons and positrons 
are selected using only the matching of a charged track reconstructed in the drift chambers and a ring pattern in the 
RICH detector. Particle identification (PID) for pions and protons is provided by the correlation between the velocity 
(/3=v/c) obtained from TOF or TOFINO scintillator walls and the momentum deduced from the track deflection 
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in the magnetic field |44j . The start signal for the time measurements was taken from the fastest signal from the 
scintillator wall. To reconstruct the time-of- flight for each particle, a dedicated method was developed [44], using the 
identification of one reference particle, the time-of-flight of which can be calculated. When a e + e~ pair was present 
in the exit channel, one of the leptons could be used as the reference particle. For events without electron candidates 
but containing a negatively charged particle, it was used as a reference particle and assigned the pion mass. When 
only two positive hadrons were observed in the final state, two hypotheses were tested: 1) the presence of two protons 
(2p events) and 2) the presence of one pion and one proton (7r + p events). For each hypothesis, both hadrons could 
be used as reference particles, hence providing an additional consistency check. In all the cases described above, the 
time-of-flight of the reference particle was calculated, and the velocities of all the other products were then deduced, 
using only the time-of-flight differences to the reference particle. The correlation between velocity and momentum of 
all particles was then used to reject the wrong hypotheses and to assign the final PID of all particles. Figure Q] displays 
such a correlation for positively charged tracks without signal in the RICH. The efficiency of the PID procedure was 
higher than 90% for both pions and protons. In addition, in the case of the four-particle exit channels, the algorithm 
was checked in a dedicated experiment with a low beam intensity using a START detector, as discussed in [4"4"] . 

2.3 Acceptance and efficiency considerations 

The spectrometer acceptance, detector efficiency and resolution as well as the analysis cuts necessary to extract 
the signal introduce important constraints on the determination of the cross sections and on the comparison of 
the experimental distributions to model predictions. For the HADES case, to compensate acceptance losses due to 
spectrometer geometry and extract 47r integrated yields, extrapolation into the unmeasured regions of phase space is 
usually achieved by means of a model. The reliability of the model to describe the shape of the relevant distributions 
also in the unmeasured regions determines the systematic errors of the acceptance corrections. 

On the other hand, direct comparisons of theoretical and experimental distributions can also be made inside the 
HADES acceptance using dedicated filters. For this purpose, the acceptances and efficiencies for the different particles 
(i.e electrons, pions and protons) were separately tabulated in matrices as a function of momentum, azimuthal and 
polar angles. The matrix coefficients have been determined using full GEANT simulations, with all reaction products 
processed through the detector, and analysed with the same programmes as done for real events. The resulting 
acceptance matrices describe the HADES fiducial volume only and can be applied as a filter to events generated by 
models. The corresponding efficiency matrices account for the detection and reconstruction process and have been used 
to correct the experimental data event-by-event. In addition, an emulator of the trigger condition was applied both 
on experimental data and simulated events. In the case of the two-hit trigger, however, the data were corrected for 
the condition of having at least one particle in the TOFINO. The detection and reconstruction efficiency was typically 
90% for protons and pions and about 50% for electrons. In addition, the yields measured in the e + e~ channels were 
corrected for the LVL2 efficiency. The latter was calculated by comparing the pp— ?>ppe + e~X yield in unbiased LVL1 
events to the yield obtained with both LVL1 and LVL2 conditions; in this way a LVL2 efficiency of 90±5% was 
obtained. 

The momentum resolution parameters were determined from the simulations in bins of momentum and polar 
and azimuthal angles and rescaled to match the resolution determined experimentally using the elastic pp scattering 
|44j . The acceptance matrices and resolution parameters, necessary to filter and smear the model generated particle 
distributions before comparing to HADES efficiency corrected data, are available from the authors upon request. 

2.4 Normalization procedure using elastic scattering 

The normalization of the experimental yield is obtained using the analysis of events produced in elastic scattering. 
Momentum conservation in the pp elastic scattering reaction leads to the two following relations between the polar 
angles 9\ and 62 and azimuthal angles cf>i and </>2 of both protons: 



where 7cm is the Lorentz factor of the center-of-mass system. The elastic events were selected by an elliptic cut in 
the (\4>i — 02 1 j tan#itan#2) plane, with semi axes corresponding to approximately 3a for each variable, i.e. ±2.4° 
for \4>i — 02 1 and a= 0.027 for tan^i tan 62 ■ The angular distributions of the resulting event ensemble is corrected for 
efficiency and compared to a simulation which uses the high precision data from the EDDA experiment [49] as input, 
as shown in fig.[5]for the 1.25 GeV incident energy case. The histogram in fig. [5] represents the angular distribution of 
events generated by PLUTO [SUllBTj according to the parameterization of EDDA data and subjected to the HADES 



\4> x = 180°, 

tan 6\ tan 62 = — — 

7cm 



(1) 

(2) 
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Fig. 2. Angular distribution in the laboratory system of measured pp elastic events (full dots) compared to PLUTO simulations 
(full histogram) using the angular distributions of [49] as input. 



filter, as described in more details in sec. 12.31 The shape of the angular distribution is well reproduced, demonstrating 
that the angular dependence of the efficiency correction is under control. The cut for angles larger than 62° reflects 
the cut on the forward partner at about 18°, which is due to the detector acceptance. The experimental yield was 
scaled in order to reproduce the simulated yield inside the HADES acceptance. The resulting factors are used for the 
normalization of the differential cross sections and have a precision of about 6% at 1.25 GeV and 11% at 2.2 GeV, 
reflecting mainly the uncertainty on the global efficiency of the reconstruction and analysis. 



3 Data analysis 
3.1 Hadronic channels 



3.1.1 Selection of pp— ^pp7r° and pp— ?>pn7r + channels 

To study the pp— >pp7r° and pp— >~pmr + channels, only events with two protons (2p) or one proton and one ir + (p7r + ) 
have been considered, respectively [5"2"1[5"3"] . The selection of both channels is based on the requirement for the missing 
mass to the system of the two detected charged particles to be close to the missing neutral particle mass. For the events 
with two detected protons, the distribution of the squared missing mass to the two-proton system (M^ liss (p,p)), shown 
in fig. [31 present for both energies a prominent asymetric peak close to zero. This contribution, clearly due to the 
elastic scattering, nicely fulfills the corresponding angular correlation (see sec. l2.4p . as shown by the dashed histogram. 
The different widths of these peaks at both energies as well as the better separation of the one-pion contribution at 
1.25 GeV result from the resolution on the proton momentum. In addition, these spectra reflect the larger phase space 
available for inelastic processes at 2.2 GeV. The contribution of the 27r contribution, visible for M^ niss (p,p) larger than 
0.08 (GeV/c 2 ) 2 , is clearly enhanced and the rj meson production shows up for M 2 Tliss (p,p) around 0.3 (GeV/c 2 ) 2 . 

To proceed with the selection of the exclusive one-pion production channels, the elastic events, selected as explained 
above, were first removed from the 2p sample. The resulting M£, isa (p,p) spectra, shown in fig. [4^ and [4]:, are peaked 
close to the squared pion mass (m 2 =0.02 GeV 2 /c 4 ). In the case of the p7r + events, the distributions of the missing 
mass to the p7r + system (M m i SS (p,7r + )) are shown in fig. Ht>, and the unmeasured neutrons become visible as peaks 
around 0.94 GeV/c 2 . The contribution of two-pion production is seen on the right hand side of the peaks. At 2.2 GeV, 
this channel was simulated as resulting from a double A production, with normalization adjusted such as to fit the 
data at the highest missing masses, as shown in fig. \%}p andHJl. For each phase space bin considered in the analysis, the 
two-pion contribution was subtracted and the remaining yield was then fitted with a function consisting of the sum of 
two gaussians plus a polynomial background. The signal was defined as the yield above this background; systematic 
errors were of the order of 5%. For the 1.25 GeV data set, the two-pion production region was simply excluded from 
the fit. Systematic errors for the signal yield were estimated from a variation of the background parameterization to 
be of the order of 1-3%. 
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Fig. 3. Squared missing mass spectra for the reaction pp— >ppX at 1.25 GeV (top) and 2.2 GeV (bottom). The dashed histogram 
shows the events selected by the elastic scattering angular correlations eqs, (IH2p . 
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Fig. 4. Selection of events from pp— >-pn7r + and pp— >pp7r° reactions based on a fit of missing masses. In the case of 2.2 GeV 
incident energy, the 27r production deduced from the simulation (dash-dotted curves) is subtracted before the fit. The fitting 
function consists of a sum of one polynomial (dotted curves) and two-Gaussians (the sum of which is shown as dashed curves). 
The solid curves show the sums of all the contributions. 
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3.1.2 Selection of the pp— !-pp7/ —> pp7r + 7r ir° channel 




M miss (p,p,7t + x) [GeV/c 2 ] 

Fig. 5. (Color on-line) Analysis of the pp7r + 7r~ events in the pp reaction at 2.2 GeV. The correlation between the missing mass 
to the two-proton system (M m i ss (p,p)) and the missing mass to the four-particle (M m ; ss (p,p,7r + ,7r~)) is shown. Events from 
the pp— »pp7r + 7r - reaction are visible for M m i SS (p,p,7r + ,7r~) close to zero. The rj peak is clearly visible for M m i ss (p,p) around 
the rj mass and M^[ sa n close to the ir° mass. The orange rectangle indicates the region used to further extract the rj signal. 

To investigate the pp— > pprj — > pp7r + 7r _ 7r° channel, events with two protons, one positive and one negative pion 
(pp7T + 7r - events) are considered [5"4lf5"5"ll5()] . Two observables have been denned: the missing masses M miss (p,p) and 
M m i SS (p,p,7r + ,7r - ) to the two-proton and four-particle systems, respectively. The correlation between these two missing 
masses is displayed in fig. [5] The concentration of events with a missing mass to the pp7r + 7r" system slightly above zero 
is due to the pp— >pp7r + 7r~ reaction; the broad structure with M m i SS (p,p,7r + ,7r _ ) around the tt° mass corresponds to the 
pp— >pp7r + 7r~7T° final state and contains an elongated spot around M miss (p,p)=0.55 GeV/c 2 , clearly due to the pp— >~ppr] 
signal. To extract the latter, first a selection of M„ l i SS (p,p,7r + ,7r _ ) between 0.05 and 0.25 GeV/c 2 , (corresponding to 
the vertical edges of the rectangle in fig. [5]), was applied in order to reject most of the pp7r + 7r~ background. 




M miss (p,p) [GeV/c 2 ] 

Fig. 6. Distribution of the missing mass to the two-proton system (M m i S3 (p,p)) for the pp7r + 7r~ events measured in the pp 
reaction at 2.2 GeV, after selection on the missing mass to the four-particle syatem (M m i ss (p,p,7r + ,7r - )) around the 7r° mass 
(full black dots). The data are normalized to the pp elastic yield. The dashed curve shows the fit of the non resonant three- pion 
background. The empty circles result from the subtraction of this background and define the rj signal. The full histogram is the 
result of the simulation of the pp— >ppr] reaction. 

The resulting (M m i SS (p,p) spectrum, normalized to the elastic yield is displayed in fig. [B] and shows a peak at the 
mass of the r\ meson on top of a broad continuum, which is mainly due to the non-resonant 7r + 7r~7r° production. Its 
contribution in the peak region (i.e. missing masses between 0.490 and 0.610 GeV/c 2 , delimited by the horizontal edges 
of the rectangle in fig. [5J was obtained from a polynomial fit of the data outside the peak region. The r\ signal was 
defined as the yield above this background, corresponding to about 24800 counts. The sensitivity to the background 
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suppression was studied by varying the limits for the fit. It gave a systematic error of the order of ±4%. The missing 
mass distribution obtained from the simulation of the pp— ^pp?7 channel is shown as a full histogram in fig. [6] Its width 
depends only marginally on the ingredients of the model for the 77 production, which will be discussed in more details 
in sec. HI The agreement of simulation and experimental signal confirms the consistency of the extracted rj signal and 
the good description of the detector resolution in the simulation. 



3.2 Exclusive dielectron channels at 2.2 GeV beam energy pp^-pp7r /^ -» ppe+e 7 




M 2 miss (p,p,e + ,e ) [GeV 2 /c 4 ] 

Fig. 7. (Color on-line) Analysis of the ppe + e~ events in the pp reaction at 2.2 GeV. Correlation between the square of the missing 
mass to the two-proton system (M^, iss (p,p)) and the square of the missing mass to the four-particle system (M^ liss (p,p,e + ,e - )) 
(see text). Events from the pp— >ppe + e~j reaction are visible for M^ liss (p,p,e + ,e~) around zero. The orange rectangles show 
the limits used to extract the 7r° and 77 Dalitz decay signals. 



x1CT 

I pp -> ppe + e"X l 

E=2.2 GeV 




Fig. 8. (Color on-line) Distribution of squared missing mass to the two-proton system M^ nis3 (p,p) (black dots) for ppe + e~ 
events after a cut on the missing mass to the four-particle system (|M^ iss (p, p, e + , e - ) | < 0.023 GeV 2 /c 4 ). The yields are 
divided by the pp elastic scattering yield. The histograms show the results of GEANT simulations. The red histograms peaking 
at the 7T° and 77 squared masses correspond to the exclusive 7r° and r\ production, respectively, followed by Dalitz decay, with 
cross sections as listed in tabled The blue histogram shows the contribution of multipion background, which is subtracted from 
the data, as explained in the text. The vertical lines depict the limits used to extract the experimental signals strength. 



To reconstruct the pp— >p-pir°/i] — !• ppe + e~7 channels at 2.2 GeV beam energy, all events with two protons and 
one dielectron (ppe + e~) have been selected, utilizing conversion pair rejection, as described in sec. 12.21 and [52ll55l 
I551I57| . The same analysis procedure was also applied to events containing like-sign pairs (ppe + e+ and ppc~c~). The 
combinatorial background, defined as the arithmetical mean of the corresponding distributions for the latter two event 
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samples, were then subtracted from the unlike-sign pair distributions. The signal-to-background ratios are of the order 
of 3 in the tt° region and of 4 in the r\ region. In the same way as for the pp7r + 7r _ events, the missing masses to the 
two-proton system M m i SS (p,p) and to the four-particle system M m j SS (p,p,e + ,e~) were reconstructed, respectively. 

The correlation between the squares of both missing masses is shown in fig. [7] after combinatorial background 
subtraction. Events from the pp— >ppe + e _ 7 reaction are visible for M^ liss (p,p,e + ,e _ ) around zero. The contributions 
from pp->pp7r°/?7 reactions followed by Dalitz decays ir°/r] — >c+c~7 can be seen for M^ iss (p,p) close to the n° and t] 
squared masses, respectively. The regions where both signals are extracted are shown as rectangles in fig. [7] The main 
remaining background is due to e + e~ pairs from the Dalitz decay of a 7r° produced in multi-pion production processes. 
Its contribution, of the order of 5% has been simulated using the cross sections 1.09 mb for 7r°7r° production and 0.50 
mb for the 7r 7r°7r°and 7r°7r + 7r _ [2], and is removed bin by bin, as illustrated in fig. [8] Systematic errors of the order 
of 3% and 8% have been estimated for the ir° and r\ signal yields, respectively, by varying the missing mass limits and 
the shape of the multipion background. 

A total amount of 6800 ± 82 stat 7r° — > je + e~ events and 235 ± 19 s t a t V — > 7e + e~ events have been extracted. 
The dashed histogram in fig. [5] is the result of the simulation of 7r° and rj Dalitz decays, with ingredients based on a 
resonance model, as will be explained in the following. The widths of the missing mass peaks, which do not depend 
on the details of the model and mainly reflect the momentum resolution of the particle tracks, are similar to the 
experimental ones. In the same way as for the other analysis channels, this gives confidence on the reliability of the 
procedure to extract the proper signal selection cuts from the simulation. 



4 Simulation of the reaction channels 

The inputs of our simulation are inspired by the resonance model by Teis et al., which is the basis for the coupled 
channel BUU transport code (CBUU) [3T]. The contribution of .4(1232), which is dominant for the pion production at 
the lowest energies, is taken from the One-Pion- Exchange (OPE) model of Dmitriev et al. [32], which describes quite 
well the measured cross sections, invariant mass distributions, and angular distributions of the pp— >pn7r + reaction at 
incident energies between 0.97 and 3.2 GeV [3"ll4"lfl"8"] . One important parameter of the model, which had been fitted 
to reproduce these data, is the cut-off parameter 71^=0.63 GeV, entering the wNA and irNN vertex form factor 

F(t) Al ~ ml 



Al-t 



, with t being the four-momentum transfer squared and m w the pion mass. The available cross section values for 
exclusive one-pion, two-pion or rj production in pp and pn reactions were used to fit the contributions of isospin 1/2 
(N*) and isospin 3/2 resonances other than the 4(1232) |31j . 

For the simulation of the channels analysed in our experiment, we employed the event generator PLUTO |50) and 
included in its data base the cross sections for the different reactions (see table [2]) . The cross sections in the two first 
columns are directly taken from [21] and the ones in the last columns are adjusted to better describe the present data, 
as will be shown in sec. The following relations derived from the isospin coefficients, are fulfilled in the simulation: 

er(pp -> 4 ++ n -> 7r + pn) = 9cr(pp -> 4 + p -> 7r + np) (3) 

g 

= 2 CT (PP A+ P n PP) ( 4 ) 

and cr(pp -» 4N — s- 7r + pn) = 5cr(pp -> 4 + p -> 7r°pp). (5) 

In the same way, one gets for the 1=1/2 resonances, 

cr(pp -> N*p -> 7r + np) = 2ct(pp -> N*p -> 7r°pp). 

Resonances heavier than N(1535), which, in the original Teis fit [31], contribute 7% and 11% to the pn7r + and pp7r° 
final states, are neglected in our approach. As described in more detail in [50], the resonance mass distributions were 
taken according to [31 . Besides the already mentioned case of the 4(1232), the angular distributions for the production 
of the other resonances are assumed to be isotropic in the pp center-of-mass frame, as in the original Teis model [3T] , 
except for the N(1440) resonance, where a steep distribution following the One-Boson-Exchange (OBE) model of [59"] 
was implemented. The decay angular distributions were kept isotropic, as in [3T], except for the 4(1232). In this 
case, the angular distribution of the 4 decay (4— >N7r) behaves as 1+0.65 cos 2 0, where is the angle between the tt 
momentum in the 4 rest frame and the momentum transfer calculated in the rest frame of the excited nucleon. Such 
a shape was indeed found to reproduce the available data [5][8]. 

For the rj production, a non-resonant contribution was introduced, in addition to the N(1535) (see tabled with 
the same proportion as in the analysis of the DISTO data 301, measured at similar beam energies, and was simulated 
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Fig. 9. (Color on-line) Dalitz plots of the pp— >pp7r° (left panels) and pp— »pn-7r + (right panels) reactions: nir + and p7r° invariant 
mass squared distributions at 1.25 GeV (upper row) and 2.2 GeV(lower row). In panel (c), the region affected by the final state 
interaction of the reaction pp— >pniv + is marked by a circle. 



following phase-space distributions. For both resonant and non-resonant contributions, the angular distribution was 
deduced from the DISTO data. The r\ production cross section, which was not measured in the DISTO experiment, is 
taken from [31] . The description of the Dalitz decay of 77 and 7T° mesons uses Vector Meson Dominance (VMD) model 
form factors (cf. [50]). 

Finally, the pp and pn final state interactions have been implemented using the Jost function to weight the 
distributions in the simulation |60j . This model, including the changes with respect to the original resonance model 
[31] mentionned above and summarized in table [3] is called "model A" in the following. In the course of the paper, a 
second model (" model B" ) will be introduced, which consists in a better parameterization of the measured data. 



5 Results and comparison with resonance model 

5.1 Dominance of A resonance in one-pion production channels 

After selecting events from the one-pion production channels, following the procedure described in sec. 13.1.11 we first 
investigate the Dalitz plots (fig. with respect to the footprints of a resonant behaviour of particle production. For 
the pp7r° channel, an accumulation of yield for Mf ntI (p,7r°)=1.5 (GeV/c 2 ) 2 , corresponding to the excitation of the A + 
resonance is clearly seen at both incident energies. 

For the pp— !-pn7r + reaction, the A ++ signal stands out markedly at M 2 nt ,(p,7r + )=1.5 (GeV/c 2 ) 2 , while the A + 
signal located at M 2 ra ^(n,7r + )=1.5 (GeV/c 2 ) 2 is not visible. The dashed curves in fig. [9] indicate the kinematical limits 
of the Dalitz plot for the different channels. The empty zones in the plots are due to the acceptance cuts, the dominant 
effect being due to the minimum proton polar detection angle of about 18°. For the pp— )-pn7T + reaction at 1.25 GeV, 
the enhanced population for both M 2 lt ,(p,7r + ) and M 2 nt) (n,7r + ) around 2 (GeV/c 2 ) 2 is due to the pn Final State 
Interaction (FSI), which enhances events with small relative momentum between the proton and the neutron. The 
FSI is less apparent in the pn7r + channel at 2.2 GeV, since it affects events with proton angles below the acceptance 
limit. For the pp7T° channel, the pp FSI has a maximum effect when both protons hit the same sector of the HADES 
detector, which is suppressed by the trigger configuration. 

Figure [10] exhibits respectively the p7r° invariant mass for the pp— >pp7r° reaction in the left part and the p7r + 
and n7r + invariant masses for the pp— >pn7T + reaction in the right part. The data are corrected for reconstruction 
efficiencies and normalized using the total pp elastic cross section, as explained in sees 12.31 and !2. 41 Error bars include 
statistical and systematic errors due to signal selection (1-5%) and efficiency corrections (5-10%). In addition, both 
isospin channels at a given energy are affected by the same global normalisation uncertainty of the order of 6% at 1.25 
GeV and 11 % at 2.2 GeV. The Mj n „(p, ir + ) and Mj n „(p, ir°) distributions are peaked around 1.23 GeV/c 2 , which 
confirms that most of the pions are produced via A decay, although the distributions are obviously distorted by the 
acceptance. The different contributions of the simulation with cross sections taken from [3T] are shown in fig. [TU] too. 
At both energies, the Mi ntI (p, tt + ) and Mi nv (p, ir°) distributions are mainly sensitive to the A contributions. The 
trend of the data is rather well reproduced, although obvious discrepancies concerning both the yields and the shapes 
can be observed. 

At 1.25 GeV, the model A overestimates the experimental yield by 20% for pn7r + and underestimates it by 20% 
for pp7r°. For the pn7r + channel, this discrepancy is slightly larger than expected by taking into account, on the one 
hand, the discrepancies of the fit to previous data in both isospin channels and, on the other hand, the combined 
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Fig. 10. (Color on-line) 7rN invariant mass distributions (full dots) measured in pp— >pp7r° and pp— >pmr + reactions at 1.25 
GeV (top row) and 2.2 GeV (bottom row). The data are compared inside the detector acceptance on an absolute scale to the 
predictions of the model A (see sec. |4| with contributions of Zi + (1232) (pink dotted curve), /A ++ (1232) (dashed blue curve), 
N(1440) (green short dash-dotted curve) and N(1520)+N(1535) (long dash-dotted light brown curve). The long dashed curve 
shows the result of the model B with a scaling factor of 0.85 applied in the case of pp— >-pn7r + at 1.25 GeV (see sec. 15. 2. 2"]) . 



uncertainty of normalisation (about ± 7 %) and global efficiciency corrections (about ± 8 %). However, the yields are 
obtained here in a limited region of the phase space, and they are therefore sensitive to the distributions used in the 
model, as will be shown in the following. 

At 2.2 GeV, the contributions of higher lying resonances clearly show up at high invariant masses and are under- 
estimated in the simulation using the cross sections from [31] (see table [2]). 



5.2 Analysis of the pp^pn7r+ channel 

We will now discuss in some detail the distributions obtained in the different channels. The results arc first compared 
to the resonance model [3T] to show its capability to describe the data and then a better parameterization of the data 
is proposed. 



5.2.1 A resonance angular distributions 

The neutron angular distributions in the center-of-mass system measured in the pp— >pn7r + channel at both energies 
are displayed in fig. 1111 These distributions mainly reflect the angular distribution of A resonance production, since 
pp— >nA ++ is the dominant process for the exclusive tt + production. They are strongly forward/backward peaked, as 
expected for the characteristic peripheral production of the A resonance. The distribution in fig. ITTb , (pp— >pn7r + at 
1.25 GeV) is highly distorted in the backward hemisphere, which is mainly due to the limited acceptance for protons 
at small laboratory angles (9 < 18°). At 2.2 GeV (fig. Illb). these acceptance losses are even larger. Therefore, we did 
not use the backward hemisphere at this energy. We included into fig. [TT] the results from the simulations. In order 
to compare the shapes of the neutron angle distributions, the simulations were rescaled to reproduce the integrated 
experimental yields. At 1.25 GeV, it can be seen that the forward/backward asymmetry is quite well reproduced by 
the simulation. The distribution is somewhat less peaked in the case of the A + and N* excitations, since the neutron 
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Fig. 11. (Color on-line) Angular distributions of the neutron in the pp center-of-mass system. Top: E=1.25 GeV, Bottom: E=2.2 
GeV. Data (black points) are compared to simulations with model A (solid curve) including Z\ ++ (blue dashed curve), A + (pink 
dotted curve), N(1440) (green dash-dotted curve), N(1520) and N(1535) (brown dot-dot-dashed curve). The dot-dashed curves 
show the result of model A with yl^=0.75 GeV in the case of 1.25 GeV and model B in the case of 2.2 GeV. Both calculations 
are scaled to reproduce the respective integrated experimental yield. 



comes from the decay of the resonance. The A ++ contribution however still dominates by an order of magnitude 
around cos 9 n =Q. The slope at forward angles is well described by the sum of the different components (solid curve), 
but the experimental distribution is slightly less steep than the simulated one, at both energies. With the chosen 
normalisation to the integrated yield, the experimental yield at cos n =0 is larger than the simulated one by factors 
of about 1.6 and 2.5 at 1.25 GeV and 2.2 GeV, respectively. 

Acceptance corrected angular distributions, which are useful to provide a result independent of the detector ge- 
ometry, can only be obtained using a model. This could be done with a good precision, at 1.25 GeV only, where the 
reaction mechanism is best under control, due to the overwhelming contribution of the A resonance. The acceptance 
correction factors are calculated for different (cos 8 n ,Mi nv (p, 7r + )) cells, chosen to optimize the precision of the correc- 
tion and defined as the ratio of events from simulation in full phase space and in geometrical HADES acceptance. In 
this way the factors depend weakly on how the invariant mass and angular distributions of the A are realized in the 
model. Remaining uncertainties come mainly from the decay angle distribution of the A resonance. Due to the limited 
acceptance of our experiment, the measured distributions of the tt + emission angle do not allow to improve the results 
from previous measurements 015], in which a decay angle distribution compatible with l+Bcos 2 # with B=0.65 ±0.30, 
was measured. The uncertainty on this anisotropy parameter has therefore been taken into account to calculate the 
systematic errors. The acceptance corrected neutron angular distribution obtained in this way is shown with statistical 
and systematic errors in fig. 1121 Once corrected for acceptance, the neutron angular distribution recovers the expected 
forward/backward symmetry. The integral of this distribution gives the cross section for the pp— >pmr + reaction, as 
will be discussed in sec. 15.41 The prediction from the resonance model is also shown, on fig. I12[ after a renormalization 
by a factor 0.85 to match the integrated yield of the experimental data. The underestimation of the experimental yield 
around cos 8 n f=s is consistent with the result obtained within the HADES acceptance. 

Since the shape of the A production angular distribution in the OPE model depends on the value of the cut-off 
parameter the sensitivity of the simulation to this parameter was studied, keeping the cross sections of the different 
contributions as in [3T]. For yl„.=0.75 GeV, instead of the standard value of 0.63 GeV, the difference between model 
and experimental data at cos 9 n =0 is reduced from 40% to about 15% (see figs. fTTk and [T2|) . It is clear, however, that 
the discrepancy of the measured angular distributions with respect to the OPE model might have a different origin 
than just a refitting of the cut-off parameter. In the region of cos 6 n =0, p meson exchange could be more important 
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Fig. 12. Angular distribution of neutron in center-of-mass system after acceptance correction for the pp— >pn7r + reaction at 
1.25 GeV. Data (black points) are compared to simulations based on model A with An-=0.63 GeV (solid curve) and the modified 
version with A n —0.75 GeV (dashed curve). Both simulation curves are normalized to reproduce the integrated experimental 
yield. 



due to the higher four-momentum transfer. In addition, non-resonant contributions might have a much flatter angular 
distribution than the A contribution. Finally, the interference between the amplitudes of the different resonances are 
neglected in our description. 



5.2.2 Modifications of the resonance model 

Considering the aforementioned deviations of the experimental distributions with respect to the resonance model of 
[3"T] , some modifications were introduced to provide a better parameterization of the data. 

At 2.2 GeV, the A resonance contributions were not changed, but the cross sections of higher lying resonances 
(N(1440), N(1520) and N(1535)) were increased and a non-resonant contribution, generated with a phase-space dis- 
tribution, was added. The new cross sections are listed in the last column of table [D A better description of the p7r + 
and 7r + n invariant mass distributions in fig. 1101 and of the neutron angular distribution in fig. [IT] in the pn7r + channel 
can be obtained, as can be seen by the long dash-dotted curves on the corresponding pictures. 

At 1.25 GeV, we have chosen to keep the cross sections of [3T], but, in order to provide a parameterization of the 
data yet more precise than the OPE model with 71^=0.75 GeV (see sec. 15.2. ip . an iterative procedure was used to 
fit the A production angular distribution such as to reproduce the measured neutron angular distribution. Due to 
the dominance of the A ++ excitation in the pp— >pn7r + reaction, this angular distribution is very close to the neutron 
angular distribution presented in fig. 1121 

Applying an overall normalisation factor of 0.85, which is consistent with the different errors, (see sec. 15. 4|) . the 
yields and shapes of the invariant mass spectra are well described by this modified resonance model, as can be seen 
by the long-dashed curve in fig. [Til] Besides, the pn FSI, introduced already in model A, while affecting only a 
very small fraction of the events, is found important to reproduce the behaviour of the distribution for the highest 
invariant masses. In this way, a new parameterization of the data is proposed at both energies, called "model B". 
The modifications with respect to model A are summarized in table [31 This parameterization will be checked for the 
pp— >ppir° channel in sec. 15.31 while in the next section, the invariant masses measured at 1.25 GeV are presented in 
more detail. 



5.2.3 Invariant mass distributions at 1.25 GeV 

To further understand the contributing mechanisms in their impact on the neutron angle cos 9 n , the p7T + and 7r + n 
invariant mass spectra were studied in the forward hemisphere in five different cos 9 n bins, as shown in figs. [T3] 
and [14J The spectra are compared to the results of the simulation with model B, which takes into account the 
experimental neutron angle distribution (see sec. 15.2.2]) . As before, the model B is scaled by a factor 0.85 to reproduce 
the experimental yield integrated over the neutron angle. The evolution of the shapes of the invariant mass spectra 
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Fig. 13. p7r + invariant mass distributions measured in pp— >pn-7r + reactions at 1.25 GeV for different bins in cos 6 n (full dots), 
compared to model B (see text) rescaled by a factor 0.85, with total (full curves), Z\ ++ (1232) (dashed curves), Z\ + (1232) (dotted 
curves), and N(1440) (dash-dotted curves) contributions. 




as a function of neutron angle in the simulation is mainly due to the detector acceptance and trigger effects on the 
dominant A ++ contribution. In particular, the structure at about 1.45 GeV for cos n < 0.8 is due to the requirement 
for the pion and proton to hit two opposite sectors. Although discrepancies of the order of 25% can be observed, 
the experimental spectra are rather well reproduced. In particular, the shape of the p7T + invariant mass distribution 
around cos 9 n =0 seems to indicate that the A contribution is still dominant in this region. The introduced flattening 
of the angular distribution of the A production both compensates the missing yield around cos n =0 in the original 
model and gives better agreement of the invariant mass spectra. These distributions definitely contain rich information 
about the pion production mechanism and should be compared to more sophisticated models including interference 
effects and non-resonant contribution. Thanks to the high statistics, these detailed distributions can indeed provide 
constraints which are complementary to the results from the bubble chamber experiments 



16 



HADES collaboration: Exclusive one-pion and one-eta production with HADES 



5" 5 
E 



pp^ppTt E=1.25GeV 



O q. 0.4 




B, total 



CD 

o 



£ 0.3 







0.2 



0.1 







-1 



-0.5 







0.5 



cosGp 



,CM 



Fig. 15. Center-of-mass proton angular distribution in the reaction pp— >ppn at 1.25 GeV compared to three versions of the 
resonance model: model A (full curve), modified version of model A with A n —0.75 GeV (long dot-dashed curve), and model 
B (dashed curve). The A + and N*(1440) contributions are depicted separately for model A by dotted and short dot-dashed 
curves, respectively. 

5.3 Analysis of the pp^pp7r° channel 

We will now discuss the results obtained in the pp7r° channel and compare them to the resonance model in its standard 
and modified versions. 



5.3.1 Invariant masses 

As mentioned in sec. 15.11 model A, based mainly on the Teis resonance model [3T], underestimates the yield in the 
pp7T° channel by 20 % at 1.25 GeV and 35 % at 2.2 GeV (see fig. ITU1) . After the inclusion of the changes in the model 
B motivated by the study of the pp— >-pn7r + channel (see sec. I5.2.2p . i.e. a slight rescaling of the cross sections of the 
different channels and the use of a phenomenological angular distribution for the A production, both the yields and 
the shapes of the 7r°p invariant mass distribution are better reproduced, as demonstrated by the long dash-dotted 
curves in fig. [TUJ 

At 1.25 GeV, the change of the A production angular distribution in the model mainly results in a global increase 
of the cross sections in the HADES acceptance by 33%, with small effect on the shape of the 7r°p invariant mass 
distribution. Note that, contrary to the pn7r + channel, no rescaling is applied to the model. In contrast to the pp— >pn7r + 
case, where the whole A production angular distribution could be measured, A production in forward or backward 
angles is suppressed by the HADES acceptance in the pp— >-pp7r° channel. This is the reason of the higher yield obtained 
for the simulation with model B, where the A production angular distribution is flatter. 

At 2.2 GeV, the change of the shape of the 7r°p invariant mass distribution is induced by the increase of the 
cross sections for the higher lying resonances (N*(1440), N*(1520), N*(1535)) and the introduction of a non-resonant 
contribution. 

To summarize, the changes motivated by the study of the analysis of the pp— >pn7r + channel also improve the 
description of the yields in the pp— >pp7r° channel, which adds consistency to the procedure. To complete the study, 
the proton angular distributions measured at 1.25 GeV are investigated in the next section. 



5.3.2 Proton angular distributions at 1.25 GeV 

Even at 1.25 GeV, where pp— >pA is the dominant process, the A + resonance cannot be unambiguously reconstructed 
due to the two protons in the exit channel. However, although both, the proton coming from the decay of the A 
resonance and the scattered one, contribute, the shape of their angular distribution inside the HADES acceptance 
is mainly sensitive to the distribution of A + production angle and depends only marginally on the decay angle in 
our simple two-step model. It is therefore interesting to check whether the distribution of the proton angle allows to 
draw conclusions on the distribution of the A angle, which are consistent with the pn7r + channel. As can be seen in 
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fig- El the acceptance is limited to intermediate proton angles (about 45° < 9^ M < 155° ) in the center-of-mass. The 
experimental distribution is however clearly much flatter than predicted by the simulation based on the resonance 
model A (see sec. 2]). Changing the A n parameter in the vertex form factor from 0.63 to 0.75 GeV, as motivated by the 
analysis of the pp— >pn7r + reaction, the simulation comes closer to the data, although the yield around cos 8 P =0 is still 
too low. This is related to the remaining underestimation of the cos#„ distribution in the pp— >pn7r + channel. A better 
agreement can indeed be obtained with the model B, which uses as an input for the A production angular distribution 
the distribution fitted to the measurement in the pp— >pn7r + channel (see sec. I5.2.TT) . as shown by the dashed curve in 
the picture. This confirms that the two isospin channels can be described consistently with the same A production 
angular distribution. Thus, model B can be exploited for the analysis of the exclusive pp— >ppe + e~ channel at 1.25 
GeV, where a realistic model for the pp— >pA + reaction is needed. 

5.4 Exclusive one-pion and one-eta production cross sections 



E [GeV] 

1 2 3 456789 10 11 

i — 1 1 i i i i i i i i i i i i i i i i i i i i i i i i 1 1 1 1 1 1 1 1 1 1 1 1 1 




\fs [GeV] 

Fig. 16. Cross sections measured by HADES (full dots) in hadronic channels for the pp— >pn7r + (top), pp— >-pp7r° (middle) and 
pp— >ppr] (bottom) reactions compared to existing data (empty dots 2 : 20 , 6lj, empty triangles [7], full triangles [H,15 ). The 
curves display the resonance model cross sections [3T] (total: full curve, Zi(1232): dashed curve, 1=1/2: dotted curve, 1=3/2 
other than Zi(1232): dashed-dotted curve). 

The cross sections for the different one- meson production channels are reported in the first row of table [4] For the 
pp— >-pn7T + channel at 1.25 GeV, the cross section was obtained by integrating over cos 9 n the acceptance corrected 
neutron angular distribution (see sec. 15.2.11) . For the other channels, the modified resonance model (model B) was 
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used to extrapolate the measured yields to 4tt. More precisely, the acceptance factors were calculated as a ratio of 
the number of events generated in the simulation to the number of events after filtering by the HADES geometrical 
acceptance and analysis cuts. To estimate the model dependence of these corrections, parameters of the model were 
varied, especially the A decay angular distribution for the pion production channels and the proportion of resonant 
contribution in the case of the r\ production channel. The corresponding numbers are indicated in the second row 
of table |3J The main sources of errors are the model dependence of the acceptance corrections (second row), the 
normalization procedure (third row) and the efficiency corrections (fourth row). The errors due to the event selection, 
following the procedures discussed in sec. I3.1.T1 and l3~TT2l are also indicated (fifth row). Statistical errors (sixth row) 
are negligible. Note that the contribution of the error on the r\ branching ratio into 3 pions is also of the order of 1%. 
The obtained cross section values are compatible with previous measurements, as can be seen from fig. 1161 

For the pion production channels, most of the data points obtained for y/s between 2.0 and 2.4 GeV come from 
KEK [7 (black dots) and were not included in the Teis fits, which were based on CERN data tables [2J. For the 
7T+ production, the KEK points [7] fit rather well with the Teis curve as well as with previous data [flfTB]. while 
they are about 15-20% higher for the ir° production. Our data are compatible with the Teis curve, despite a slight 
underestimate of the 7r° production at 2.2 GeV. As already mentioned, the relatively large error bars of our data 
with respect to the existing previous data are due to the combined effects of efficiency corrections, normalization and 
acceptance corrections, which were reduced in the case of bubble chamber experiments. The scattering of this data 
collection might however point to an underestimate of the error related to the event identification. From the HADES 
measurements, the ratios of pp— >pn7r + to pp— >pp7r° cross sections 4.57±0.54 at 1.25 GeV and 3.49±0.63 at 2.2 GeV 
can be deduced, which has to be compared with the factor 5 expected in the case of A excitation only, see eq.©. 

As for the 77 production, our experiment brings a new measurement (a = 0.142 ± 0.022 mb) of the exclusive 
production cross section in pp reaction, in a region, about 230 MeV above the threshold, where only the two mea- 
surements from Pickup et al. |62j existed. Our point is in agreement with their value obtained in neutral channels, 
a = 0.197 ± 0.077 mb, while it is found below their more precise measurement obtained in the three-pion channel 
[a = 0.242 ± 0.043 mb). The quoted error might however be underestimated, as discussed in [3D], considering the un- 
certainty due to the non resonant background subtraction. The cross section parameterizations used in [30 and based 
on fits of data with yfs ranging from threshold up to 3.4 GeV did not take into account these points and provided, 
at an energy of 2.2 GeV, values between 70 and 100 fib, which are much lower than both Pickup's results [62] and 
the value measured in the present experiment. Our measurement is in very good agreement with the resonance model 
[3T] . where the r/'s are assumed to be produced only via N(1535) resonance decay. This assumption of fully resonant 
production seems however in contradiction with the DISTO analysis [30] . Our new measurement can hence be used to 
test various models of r\ production. Previous OBE calculations [63,64,65 showed deviations of a factor 2 depending 
on the values of the pNN*(1535) and wNN*(1535) coupling constants, which should be updated in view of the actual 
constraints on these parameters, as was done recently closer to threshold [5B] . 



5.5 7T° and 77 Dalitz decay analysis 

5.5.1 Dielectron invariant mass 

As already explained in sec. 13.21 the 7r° and r/ Dalitz decay signals has been extracted, in each e + e _ invariant mass 
bin, using the missing masses to the pp and ppe + e _ systems and were then efficiency corrected. In addition, an 
acceptance correction obtained from simulations with the resonance model was applied. The resulting c + e~ invariant 
mass distributions are displayed in fig. 1171 with statistical errors and systematic errors added quadratically. In the 
case of the 7r°, the largest source of systematic error is the rejection of e + e~ pairs from photon conversion, while in 
the r\ region, it is due to e + e~ pairs from 7r° decay in multipion production processes. 

The experimental values are compared to the results of the simulation, with exclusive meson production cross 
sections from table [2] and branching ratios from table [1] The good agreement obtained for both the ir° and ij peak 
is therefore a check of the consistent extraction of the corresponding signals, which is very useful for all dielectron 
analyses performed with the HADES detector. The small excess around 0.03 GeV/c 2 is most likely due to a remaining 
contamination of conversion pairs. The possible contribution of Dalitz decays of baryon resonances, corresponding 
to a ppe + e~ final state, has also been investigated and is found to be negligible, except in the mass region close to 
the kinematical limit (Mi nv (e + ,e~ ) = 0.547 GeV/c 2 ), which could possibly explain that the measured yield for the 
reaction pp— is higher than the simulation above 0.5 GeV/c 2 . 

It has been checked that the shapes of these spectra do not depend on the ingredients of the simulation related to 
the meson production mechanisms, like the relative yields of the different resonant contributions, but are characteristic 
of their Dalitz decay. The description of these Dalitz decay processes in the simulation implies electromagnetic form 
factors which can be implemented in the simulation following the VMD model |501I67|. Modelling the transitions as 
point-like (refered to as QED) or using VMD form factors lead to negligible differences for the ir° — > -fe + e~ case and 
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Fig. 17. Distribution of di-electron invariant masses for the it (left) and r\ (right) Dalitz decays, obtained after efficiency and 
acceptance corrections and compared to the simulations with model A (full line), including VMD models. The yields have been 
divided by the elastic scattering yields. In the case of the 77, the simulation without r\ form factor (labeled as QED) is shown as 
a dashed curve for comparison. 



show up for the r\ — > je + e~ case only at larger values of the e + e~ invariant mass. There, our data are however not 
precise enough to provide any further quantitative constraint to these models. 

The yields are well reproduced by the simulation with meson production cross sections from table [2] In the case of 
the 77 production, since the cross section is only fixed by our measurement in the hadronic channel (sec. 15.41) . this shows 
the consistency of the hadronic and leptonic reconstructions and the good control of the corresponding efficiencies. 
More quantitatively, the ratio of yields measured in 7r + 7r~7r° and 7e + e~ decays of the 77 meson is 218±25, i.e. fully 
consistent with the value 230±5 given by the simulation. In the case of the ir°, where the production cross sections 
are constrained by independent data, the analysis of the 7r° — > ^ye + e~ channel provides a global consistency check of 
the whole analysis chain for dileptons. 



5.5.2 Helicity angle 

An interesting feature of the Dalitz decay of pseudo-scalar mesons is the polarization of the virtual photon which is 
transverse. As a consequence, the distribution of the helicity angle a e7 * follows l+cos 2 a e7 » . The calculation of this 
angle first implies a boost of all particles in the meson rest frame. Then the helicity angle is defined as the polar angle 
of the electrons in the virtual photon rest frame, with respect to the virtual photon direction. The acceptance and 
efficiency corrections were calculated using the simulation of the dielectron production via rj Dalitz decay, as described 
above. As shown in fig. 1181 this angular distribution can be fitted by a function of the form a(l+bcos 2 a e7 * ) with 
b = 0.98 ± 0.48, in agreement with the QED prediction, b=l |45j . 

With the HADES set-up, it is therefore possible to reconstruct the helicity angle distribution of the 77 Dalitz decays. 
The extraction of the anisotropy parameter b should then be also possible in the case of the A Dalitz decay, where, 
the polarization of the virtual photon is also mainly transverse, since the Coulomb amplitude in the NA transition 
is small and hence a l+cos 2 a e7 * distribution is expected. Helicity angle distributions have also been investigated in 
heavy-ion reactions [32J in order to identify the nature of the "excess" beyond the 77 contribution. 



6 Summary and outlook 

HADES has provided a measurement of the reactions pp— >-pn7r + and pp— >-pp7r° at 1.25 GeV and 2.2 GeV and pp— >pprj 
at 2.2 GeV using both hadronic and leptonic channels. Using the hadronic channels, high statistics differential cross 
sections could be measured in the HADES acceptance. In addition, integrated cross sections were extracted for all these 
channels and the neutron angular distribution in the pp— ?>pn7r + reaction at 1.25 GeV was fully reconstructed. These 
data allow to test pion production mechanisms and the contribution of baryonic resonances with a high statistical 
precision, in complement to previous low-statistics but high-acceptance experiments. We left for further studies the 
comparison of these data to calculations including resonant and non resonant contributions in a coherent way. Our aim 
in this paper was twofold: first, to show the sensitivity of the present data to the ingredients of the transport models 
used for the dielectron production, which are based on resonance models and, second, to obtain a parameterization of 
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Fig. 18. Angular distributions of electrons and positrons in the virtual photon rest frame (helicity) after acceptance and 
efficiency correction, using an arbitrary normalization. The solid curve is the fit using (a(l+bcos 2 a e7 * )) with b=0.98±0.48. 



meson and baryon resonance production for dielectron channels analysis. Following this line, an analysis based on a 
resonance model [31] was presented. An overall agreement with the original model is shown, but a better description 
could be obtained in both isospin channels, using at 1.25 GeV a less steep angular distribution for the 4(1232) 
resonance production and at 2.2 GeV an increased production cross section for the higher lying resonances. A precise 
description of the 4(1232) production angular distribution at 1.25 GeV is especially important for the on-going analysis 
of the Dalitz decay of the 4(1232) resonance using the pp— >ppe + e~channel. On the other hand, further information on 
higher lying resonances can be gained by studying two-pion production channels, which were also recently measured in 
the HADES experiments. The present determination of the exclusive 77 production cross section is most important, as 
it provides the first precise measurement of the exclusive production cross section in a region where deviating model 
predictions can be found. 

The reconstruction of 7T° and r\ Dalitz decay signals presented in this paper is fully consistent with the hadronic 
channels, and the invariant masses and acceptance corrected helicity angle distributions are in good agreement with 
QED predictions. These results confirm the ability of HADES to reconstruct sensitive observables in dielectron chan- 
nels, which is a very important consistency check for previous and next-coming analyses. The helicity angle was used to 
study dielectron sources in heavy-ion reactions |42) in different invariant mass regions and is also used to discriminate 
the 4(1232) Dalitz decay process from the pp Bremsstrahlung contribution in the on-going analysis of the exclusive 
ppe + e~ channel in pp reactions at 1.25 GeV [68] . 

As a final conclusion, the present analysis provided important consistency checks for dielectron studies, as well 
as precise results for meson production measured in hadronic channels, paving the way for further theoretical or 
experimental studies of exclusive dielectron and hadronic channels in elementary reactions. 
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Table 2. Cross sections used in the simulation. Elastic pp cross sections taken from [49] and [58] are used for the normalisation 
of the measurements. For the inelastic channels, the first set of cross sections (oTeis) is taken from [3T] and is used in model 
A (see sec[4]) at both energies and in model B (see sec !5.2.2]l at 1.25 GeV. For the r\ production, the ratio of N*(1535) to 
non-resonant production is taken from DISTO [3D]. The second set (cr a( jj), used in model B at 2.2 GeV, is adjusted to the 
HADES data (see text). 





model A 


model B 


1.25 GeV 


(1) 


(1) and (2) 


2.2 GeV 


(1) 


(1) and (3) 



Table 3. Summary of the modifications introduced in models A and B with respect to the resonance model [3T] (see text 
for more details). (1): pp and pn Final State Interaction, anisotropic .4(1232) decay angular distribution, N(1440) production 
angular distribution from [59]. (2): A* cut-off parameter changed from 0.63 GeV to 0.75 GeV and 4(1232) production angular 
distribution further adjusted to describe the neutron angular distribution in the pp— >pn7r + channel. (3) change of production 
cross sections for N(1440), N(1520) and N(1535) resonances and introduction of a non-resonant contribution following table[2] 



reaction 


pp— >pn7r + 


pp^ppir 1 ' 


PP^PP'f 


energy 


1.25 GeV 


2.2 GeV 


1.25 GeV 


2.2 GeV 


2.2 GeV 


cross section (mb) 


17.1 ± 2.0 


14.45 ± 3.2 


3.74 ± 0.48 


4.15 ± 0.85 


0.142 ± 0.022 


acceptance corrections 


± 1.0 


±1.1 


± 0.2 


±0.2 


± 0.006 


normalization 


± 1.1 


± 1.6 


± 0.25 


± 0.46 


± 0.016 


efficiency 


± 1.3 


±2.5 


±0.33 


± 0.65 


±0.013 


event selection 


± 0.3 


± 0.7 


± 0.12 


± 0.2 


± 0.005 


statistics 


± 0.01 


± 0.01 


± 0.003 


± 0.004 


± 0.002 



Table 4. Cross sections for exclusive meson production channels measured by HADES in hadronic channels are given with the 
total error, calculated as the quadratic sum of the statistic and systematic errors listed in the following rows. 



